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This article presents an analytical solution of the effective index of the 
fundamental waveguide mode of ID metallo-dielectric grating for Transverse 
Magnetic (TM) polarization. In contrast to the existing numerical solution 
involving transcendental equation, it is shown that the square of the effective 
index (n^//) of the fundamental waveguide mode of ID grating is inversely 
proportional to the slit width (w) and the refractive index (n m ) of the ridge 
material and varies linearly with the incident wavelength (A). Further, it has 
also been demonstrated that the dependence of riEff on the grating period 
(P) and the incidence angle (0) is minimal. Agreement between the results 
obtained using the solution presented in this article and published data is 
excellent. © 2015 Optical Society of America 

OCIS codes: 050.0050, 260.1960, 260.2110, 260.3910. 

1. Introduction 

With the advances in nano and micro fabrication technologies, sub-wavelength structures 
are now readily achievable DEI- Due to their high brightness in resonance, recently ID 
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metallo-dielectric grating structures with sub-wavelength slits (w < A, see Fig. [T]) have been 
proposed as useful in flat panel displays, Scanning Near-held Optical Microscopy (SNOM), 
opto-electronic devices, photo-lithography and tunable optical filter pa- 

High brightness or resonance in ID grating can be explained using two different theories [5]. 
In the regime A ~ P, where A is the wavelength of the incident light and P is the period 
of the grating, coupling between surface plasmon polariton (SPP) of opposite faces of ID 
grating is responsible for enhanced transmission [5], whereas for thick enough grating and 
A » P, resonance coupling between a diffraction order and a waveguide mode plays a major 
role in the extraordinary transmission through ID metallo-dielectric grating structures ra¬ 
in the latter case, depending upon the incident wavelength, slit width w and period P, 
different waveguide modes, both propagating and evanescent, are excited inside the slits j,7]. 
Propagating modes transfer incident energy from one side of the grating to the other and 
redistribute the transferred energy among the diffraction orders |8j. As the slit width w 
decreases, more and more modes become evanescent and a very few propagating modes 
survive [§]. In particular, when w becomes smaller than A/(2 rid), where rid is the refractive 
index of the slit/groove, only the fundamental mode propagates and most of the transmitted 
energy is carried out by this mode [T1I21J5U9] . Considering this phenomenon Lalanne et al. have 
developed an analytical model of transmission through ID grating for TM polarization [T]. 
This model can accurately predict resonance wavelengths and their diffraction efficiencies [Id]. 
Finding transmission efficiency using this model requires effective index of the fundamental 
mode, which is defined as riEff = k z /ko, where k z and ko are the ^-component of the 
wave vector of the fundamental waveguide mode and the wave number of the free space 
incident electromagnetic illumination respectively [2]. This model also depends on the grating 
parameters i.e. w, P and h, where h is the thickness of the gratings. Similarly, Porto et al. [5] 
and Garcia-Vidal et al. |10j have developed models of transmission through ID gratings by 
considering only the fundamental mode and their results agree closely with those of Lalanne 
et al. [lj. Profile of the fundamental mode and x— and components of its wave vector 


2 


{k x and k z , see Fig. Q]) as well as those of other modes can be determined numerically by 
solving transcendental modal equation proposed by Sheng et al. in 1982 [7|. This method is 
known as modal analysis and its solutions, also known as eigenmodes, correspond to various 
modes of the waveguide structure. Overall this method provides exact description of the 
modes mm and is becoming popular [IMS] due to its phenomenological interpretation 
of the wave propagation via grating structure. Despite this, the method is still a numerical 
technique and, as inherent to numerical techniques, is devoid of physical insights i.e. can 
not provide a direct relationship (such as n 2 E ^ varies inversely with w) among n E ff, w, P, 
A, 6 (incidence angle) and n m (refractive index of the ridge metal) and hence the physics 
behind the wave propagation via ID gratings is not well understood. Also, finding a solution 
of the transcendental equation requires searching inside the variable domains [TS] and is 

computationally demanding 1113 - 

In contrast to the existing numerical solution, in this article we attempt to provide an ex¬ 
plicit relation involving n E ff, w, P, n m , 9 and A. This direct relationship between n E ff and 
the grating parameters provides a vivid explanation of the physics behind the wave propa¬ 
gation via ID grating structures. To the best our knowledge, this kind of analytical model 
relating n E ff and the grating parameters is nonexistent in the literature for ID waveguide 
structures even though such a relationship exists for 2D waveguides [T6] . Further, as with 
analytical solutions, finding riEff using the method presented here is easy and does not re¬ 
quire searching inside the variable domains and consequently it is computationally much less 
demanding. Also, results obtained using our model agree very closely to those of the exact 
numerical calculation. In the process of deriving our main result, we assume w < A/(2n^). 

2. Analytical model 

Let us consider a TM polarized electromagnetic wave H y = exp (ifc 0 (sin Ox — cos 0z)) * 
exp (—iut) is incident upon the metallo-dielectric grating of Fig. Q] at an incidence angle 
6. This incident wave excites various waveguide modes which in turn transfer energy from 
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Fig. 1. ID Lamellar Grating 

the incident side (z > 0, see Fig. [Q) of the grating structure to the outgoing side (z < —h ). 
x— and components of wave vectors corresponding to different waveguide modes can be 
obtained by solving transcendental Eq. © 0 . 


cos(k 0 P sin 9) — cos (/3rP) cos (afP) + - H——] sin (/3rP) sin (ctfP) = 0 (1) 

2 p e m a 

where a = k 0 ^Jed — n-Eff anc l P = ko\] e m — n Eff are the x— components of a waveguide 
mode in the slit/groove and ridge material respectively. r = (P — w)/P,f = l — r = w/P, 
6d = n'd is the dielectric constant of the slit and e m = n 2 m is the dielectric constant of the 
grating ridge. For metallic ridges, dielectric constant is given by e m = = (rj + iht) 2 , where 

r) and k are the real and imaginary components of the refractive index. When \e m \ » \n 2 E j^\, 
P can be written as P — kon m . Assuming = vP d — 1, a can be written as a = k 0 p, where 
p = 1 — n 2 E ff. Considering above Eq. ([T]) can be rewritten as Eq. (J2]) . 
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cos(/c 0 -P sill 9) - cos (k 0 n m rP) cos (k 0 pfP) + + -—] sin (k 0 n m rP) sm(k 0 pfP ) = 0 (2) 

2 P £mOL 


Given that |e m | >> |n|^| and |e m p 2 | >> 1, the first factor of the 3 rd term of Eq. ([2]) can 
be written as ~ y/^mP/2 = n m p/2. Based upon this Eq. (J2]) can be written 

as Eq. (EJ). 


cos ( k 0 P sin 9) — cos q cos p + -n m p sin q ship = 0 


(3) 


where p = 2npw/\ and q = 2irn m (P — w)/ A. Expanding ship and cosp into Taylor series, 
Eq. (ED can be expressed as Eq. (J4J. 


1 t P ,P \ n P . P P . , . /, p ■ m n 

-n m p{p- — + — - ...) sm q — + + •••) cosg + cos(fc o Psin0) = 0 


Tim 

4tt/7 


4 6 

/ 2 V , P 
(P “ 3! + 5! 


2 4 6 

sing — (1 — + ...) cos q + cos(k 0 P sin 6) = 0(4) 

2! 4! 6! 


In the typical operating conditions where only the fundamental mode survives such as 
(r = 0.8571, / = 0.1429, \n m \ « 18.73) 0 , (0.90 < r < 0.9889, 0.0111 < / < 0.10, 
\n m \ « 6.796) g] and (0.50 < r < 0.95, 0.05 < / < 0.50, \n m \ « 5.03) [IS], |g| becomes 
greater than 1 and |p| is less than unity. In this case any power of p above 4 in Eq. (ED can 
be neglected. After some simple algebraic manipulations, one can write Eq. (]4j) as Eq. Q. 


Ap 4 - Bp 2 + C = 0 


(5) 


where A = 


f Tim sin q I cos q \ p / n m sing 

V 71-/7 4! "T" 4! h 1 271-/7 2! 


+ ^),G=(cosg 


cos(koP sin 9 )). Eq. (J5D can 
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be easily solved using the standard algebra and the solutions are given in Eq. (J6]). 


V = 


B ± VB 2 - 4 AC 


n 


Eff 


= 1 - 


2 A 

3A 2 1 

47 t 2 w 2 1 + D ± 


f 2 cosg — 8 cos (k 0 P sin6 , ) 1 1 1 1 

' 3 cos q '~D + 


( 6 ) 


where D — 1 + sing and cosg can be expanded as sin {k 0 (P — w)r]} * 

cosh {k 0 (P — w)k} + i cos {k 0 (P — w)r]} sinh {k 0 (P — w)k} and cos {k 0 (P — w)r]}* 
cosh {k 0 (P — w)k} —i sin {k 0 (P — tc)r 7 }*sinh {k 0 (P — w)k} respectively. For k 0 (P—w)n > 1, 
cosh {k 0 (P — w)k} ~ sinh {k 0 (P — w)k}. In this case, sing and cosg can be written as - 
cosh {k 0 (P — w)k} * exp [i{7r/2 — k 0 (P — w)r]}] and cosh {k 0 (P — w)k} exp [ik 0 (P — w)rj\ 
respectively. After some simple manipulation one can End Eq. (JTj). 


n 


Eff 


= 1 - 


3A 2 


47 t 2 w 2 


[1 + 


7 TW 


7 tw + i\n r 


± [1 - 
1 3 


8 cos (koP sin 6) exp {iko(P — ui)r]} 


3 cosh {k 0 (P — w)k} 

7 TW 


TT 2 W 2 


7 tw + i\n m (ttw + iXn m ) 


r] 1/2 ] (7) 


Given that cosh {ko(P — w)k} » 1 and | I < 1) the term under the square root in 

Eq. ((TJ) can be expanded into binomial series. Neglecting terms of the expansion with power 
two or more, Eq. (JTJ) can be expressed as Eq. (JHJ). 


n 


Eff 


= 1 


3A 2 


:[1 + 


7 TW 


47 t 2 w 2 ttw + i\n r 


±[l-[ 


1 4 cos (koP sin 6) exp {iko(P — w)r]} 1 

3 


3 cosh {k 0 (P — w)k} 

71W 


n 2 w 2 


7 tw + i\n m 2(7 tw + iXn m ) 2 


]] ( 8 ) 


Provided that Re(n E ff ) > 1, the solution corresponding to the fundamental mode can be 
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written as Eq. (J9]). 


2 _ 1 _ [1 _ cos (koPs'mQ) exp {iko(P-w)rj} i_A^_._ 3A 2 

Eff cosh {ko(P— w)k} -I Trw^Trw+iXrim) 8{jrw-\-i\n rn )‘ 2 


= l + i[l 


cos (fco-Psin 9) exp {iko(P—w)r)} i A 

cosh {ko(P—w)K,} -I 7rnmw(l—i ) 


8 ^( 1 -*^) 


( 9 ) 


Considering 0 < 14^ I < 1 and expanding the denominators of the 2 nd and 3 rd terms of Eq. 

AYlm 

(jnj) into binomial series of and keeping only the first and second power of w and n m of 
the expansion respectively, one can represent Eq. (J9]) as Eq. (HDD. 


n E ff - 1 


[— 

8 < 


A 


+ 


A cos (koP sin 6) 


nwn m L n^ TTwrim cosh {k 0 (P — w)k] 


exp{ffc 0 (-P — w)rj} (10) 


3. Results and discussions 

The real component of nEff corresponding to the fundamental mode obtained from Eq. (fTUj) 
is plotted in Fig. [2] for = 1, 6 = 0°, P = 900 nm and A = 1433 nm for silver grating as 
a function of slit width w. Grating parameters used in this example have been taken from 
Astilcan et al. [2], where the authors show how the effective index of the fundamental mode 
evolves as the slit width of the grating changes. For the purpose of comparison Re{riEff) from 
Ref. [2] has also been included in Fig. [2j One can see that there is an excellent agreement 
between our results and those from Ref. [2]. It is also evident that as w increases, Re(nEff) 
decreases or there is an inverse relationship between Re(riEff ) and w. To confirm this let us 
look more closely at Eq. (HO- Upon consideration one can find that the contribution from 
the 3 rd term in Eq. (ITU toward n 2 E jj and hence toward Re{riEff) is very negligible since 
cosh {k 0 (P — w)k} » 1, ] exp{ik 0 (P — w)rj}\ < 1 and | cos (fco-Psin#)! < 1. Considering this 
one can rewrite Eq. (ITU as Eq. (TlTh from which it is understandable that for fixed P, n m , A 
and 6, n 2 E ff and therefore nEff vary inversely with w. At this point we quickly mention that 
this kind of physical insight is not understandable from the existing numerical solution. 
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Fig. 2. Real part of the effective index of the fundamental mode as a function 
of slit width corresponding to rid — 1, 6 — 0°, P = 900 nm, A = 1433 nm for 
silver gratings 
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From Eq. CD, one can also find that as w —* 0, n E ^ approaches infinity as the grating 
becomes impermeable to light. On the other hand, when w —> A/2, n 2 E ^ approaches a 

constant value of 1 + <5, where 5 = i —--t- 4- depends on the dielectric constant of the 

ridge material and should be much much smaller than unity as |n^J >> 1 has been assumed. 
When w = A/2 and |n^J >> 1 (which is true for most of the metals in the infrared region 
of the electromagnetic spectrum), 5 —> 0 and n 2 E ^ (accordingly n^//) approaches unity as 
expected. However, when 0 < w < A/2 and |n^| >> 1, Eq. (ITT]) can be further simplified to 
Eq. ([I2D from which one can observe that ri 2 E ^ is inversely proportional to w and n m and 
varies linearly with A. For the purpose of demonstration we have plotted data corresponding 
to Eq. (ROD and Eq. (R2D in Fig. [3] for the same set of grating parameters of Ref. [2], It can 











Fig. 3. Real part of the effective index of the fundamental mode as a function 
of slit width corresponding to rid — 1, 6 — 0°) P — 900 nm, A = 1433 nm for 
silver gratings 

be seen from Fig. [3] that there is no difference between the curves corresponding to Eq. (TTOT) 
and Eq. (fT2]h 


n 


Eff 


— 1 + i- 


TTWrir 


( 12 ) 


Moreover to verify that n 2 E ^ varies inversely with n m , in Fig. 0 we have plotted both 
Re(riEff) and Im(riEff) as a function of the refractive index of the ridge material using 
the exact transcendental Eq. © and Eq. (fT2]l . This is equivalent to considering various 
ridge materials while keeping geometrical grating parameters i.e. w, P and h intact. Grating 
parameters have been taken from HZ], where designing ID grating for extraordinary optical 
transmission is considered using the numerical optimization technique. From Fig. |4l one 
can see that there is an excellent agreement between the exact method and the simple 
analytical method we have presented above. It is also noticeable that as Im(n m ) increases, 
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Fig. 4. (a) real and (b) imaginary components of the effective index of the 
fundamental mode as a function of Im(n m ) while keeping Re(n m ) constant 
corresponding to rid = 1, 0 = 0°, w = 21 nm, P = 150 nm, A = 1500 nm. 
Grating parameters have been taken from Ref. S3 
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Fig. 5. (a) real and (b) imaginary components of the effective index of the 
fundamental mode as a function of A where rid — 1, 0 = 0°, w — 21 nm and 
P = 150 nm. n m = 0.530 + 9.50701 is that of gold at A = 1500 nm. Grating 
parameters have been taken from Ref. HU 

loss associated with the fundamental mode ( Im(riEff )) decreases. This is due to fact that 
as the imaginary component of the dielectric constant increases, a metal becomes highly 
reflective and waves can propagate without incurring much loss. Further, although A and n m 
are related, for completeness in Fig. [ 5 ] we have plotted both Re{riEff ) and Im(riEff ) as a 
function of A while keeping n m , w, P and h constant for 6 = 0°. The grating geometrical 
parameters have been taken from Ref. tea like before and n m = 0.530 + 9.5070z is that of gold 
at 1500 nm. Data for this graph have been obtained from the exact numerical calculation 
(Eq. (Op)) and Eq. (fl2]h It can be observed that Re(riEff ) and Im{riEff ) vary linearly with A 
as predicted by Eq. (fl2j) and the agreement between the exact calculation and the simplistic 
model of riEff presented in Eq. (1T2]) is very good. 

To complete the investigation of the dependence of riEff on the grating parameters, let 
us consider the impact of the two remaining parameters, namely the grating period and the 
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incidence angle on n E ff. According to Eq. (fT2j) . n 2 E ^ and consequently n E ff do not depend 
upon these two parameters. However, if one considers Eq. (HO, then it is found that the 
dependence of n 2 E jf on P and 6 is very weak. To confirm this we have plotted n E ff as a 
function of P and 6 in Fig. [6] (a) and (b) respectively. In both cases there is a very good 
qualitative agreement between the exact result obtained by numerically solving Eq. (JT|) and 
the calculation performed using Eq. (HUD . It can be observed that the dependence of n E ff 
on P and 6 is very minimal. In particular, the difference between Re(n E ff) corresponding 
to the two extreme incidence angles (9 = 0° and 6 = 90°) is approximately 0.076%. This is 
due to the fact that the denominator of the 3 rd term of Eq. (fTUD is much larger than the 
6 dependent cos (fcoPsin#) (varies between —1 and +1) in the numerator. Consequently, a 
relatively small variation in cos (koPsinO) caused by the variation in 6 does not resonate a 
significant change in n E ff. On the other hand, the dependence of n E ff on P is discernable 
up to a certain value of the grating period, beyond that it becomes independent of P. This 
behavior of n E ff in regards to P can be explained by considering Eq. (fTUD again. For a fixed 
w, n m and 6, cosh {ko(P — w)k} is a real number and is greater than 1. As P and hence 
(P — w) increases, the denominator of Eq. (fTUD gets bigger and bigger. When the value of 
cosh {&o(P — w)k} is relatively small, the contribution from the 3 rd of Eq. (fTUD towards n 2 E ^ 
is appreciable but when it becomes immensely large, the 3 rd term from Eq. (fTUD can be 
completely ignored and n E jf becomes independent of P and approaches the value predicted 
by Eq. (1T2]) . In general, it can be concluded that if 27r(P — w)n > 10A then n E ff does 
not depend on P significantly. Further, considering the above discussion, one can conclude 
that the analytical model presented in Eq. (1T2]) for the fundamental mode of ID grating 
structure is a very good representation of the exact solution and captures all the physics 
of the fundamental mode propagation via ID grating structure. Another interesting point 
that can observed from Eq. (fT2|) is that, unlike the 2D structures such as the rectangular 
waveguides [16j, ID waveguide has no cutoff wavelength above which all modes including 
the fundamental waveguide mode are non-propagating. 
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Fig. 6. Re(jiEff ) as a function of (a) P and (b) incidence angle 6 corresponding 
to rid — 1, w — 21 nm and A = 1500 nm. For (b) P is equal to 150 nm. 


Additionally, to show that the analytical solution of ue/j presented in this article is suitable 
for different ridge materials, incidence angles and geometrical grating parameters, we present 
more data from the literature and compare them with the data generated using the model 
presented in Eq. (TT2lh In this regard, we consider data from Ref. [15] where the authors 
consider aluminum gratings, from Ref. [2] where silver grating is investigated and from 
Ref. [H] in which loss less metals are considered. Comparative results are presented in 
Table CD In all cases, irrespective of ridge materials and incidence angles, the agreement 
between the current results and those from the literature is good. We present one more 
example from Ref. [6] where the authors show the negative roles of SPP on EOT for the case 
of ID transmission grating for the TM case. In their analysis the authors find the optical 
transmission of the zeroth diffraction order via ID grating by using a one-mode (fundamental 
mode) model of optical transmission [Tj. In this model effective index of the fundamental 
mode is necessary and is determined by the method of line [Tj . As per the authors analysis, 
three different transmission peaks appear at 3.58 pm, 4.9 pm and 9.5 pm corresponding to 
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a gold grating with w = 0.50 /im, h = 4.00 /im, P = 3.50 /an, = 1 and 6 = 0°. For the 
purpose of comparison we have plotted the zeroth order transmission efficiency using the 
same set up of Ref. [6] except «e// which we have determined using the solution developed 
in this article (Eq. (fT2]R From Fig.|7]it can be seen that there are three extraordinary optical 
transmission peaks at 3.57 /im, 4.9 /in i and 9.56 /im. Upon comparison with the data from 
Ref. [6], one can find that the agreement between the data generated using our analytical 
solution and those by calculating riEff numerically is very good. 


Table 1. Effective index corresponding to the fundamental mode with rid — 1- 
A, w and P are given in nm 


Grating Parameters 

n Eff 

Literature 

n Eff 

Present Work 

Ridge- Silver 

A = 1183, w = 90 

p = 900, 9 = 0° 

1.224 + 0.002i [2] 

1.220 + 0.002* 

Ridge- Unknown 

A = 632.8, w = 93.52 

p = 500, 0 = 30° 

1.105 [ig 

1.103 

Ridge- Aluminum 

A = 450, w = 100 

P = 200, 9 = 35° 

1.142 + 0.015* pR] 

1.133 + 0.013* 


Finally, we stress that the accuracy of the solution of the fundamental mode presented 
above significantly depends on the validity of the assumption \n 2 E ^\ « |e| and whenever 
this condition is not satisfied there will be a mismatch between the n E f/s calculated by Eq. 
(JIT) and Eq. (TT2T) . It is also important to mention that retaining higher powers of p in Eq. 
([5]) does not significantly improve the accuracy of riEff but increases processing difficulties 
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Fig. 7. Zeroth order transmittance corresponding to P = 3.50 //rn, w = 0.50 
/irn, h = 4.00 /irn, = 1 and 0 = 0° [6|. n m is that of gold [18]. Transmission 
efficiency is based on the model of Ref. HIE! where riEff is needed to complete 
the calculation. In Ref. [6 ]riEff has been found using a technique called method 
of line while in plotting this graph we have used Eq. (fl2|i 
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and overall n#// becomes an obscure function of w, P, 9, A and n m . Lastly, even though 
we have not considered other grating modes explicitly, conclusions similar to those of the 
fundamental mode may be applicable to them. 

4. Conclusion 

In conclusion, we have provided a simple analytical solution of the effective index of the fun¬ 
damental waveguide mode of ID grating structure for TM polarization. It has been shown 
that the square of the effective index of the fundamental waveguide mode is inversely propor¬ 
tional to the slit width and the refractive index of the ridge material. Dependence of riEff 
on the grating period and the incidence angle is negligible. The solution provided in this 
work is very easy to compute and produces results that match closely to those of the exact 
method. We have also demonstrated that irrespective of the grating materials, incidence 
angles and incidence wavelength, the analytical solution presented in this article provides 
reliable results. 
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